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ABSTRACT 

We have measured solar-like oscillations in red giants using time-series photometry from the first 34 
days of science operations of the Kepler Mission. The light curves, obtained with 30-minutc sampling, 
reveal clear oscillations in a large sample of G and K giants, extending in luminosity from the red clump 
down to the bottom of the giant branch. We confirm a strong correlation between the large separation 
of the oscillations (Af) and the frequency of maximum power (^ max )- We focus on a sample of 50 
low-luminosity stars (f max > 100 /iHz, L < 30 £©) having high signal-to-noise ratios and showing the 
unambiguous signature of solar-like oscillations. These are H-shell-burning stars, whose oscillations 
should be valuable for testing models of stellar evolution and for constraining the star-formation rate 
in the local disk. We use a new technique to compare stars on a single echelle diagram by scaling their 
frequencies and find well-defined ridges corresponding to radial and non-radial oscillations, including 
clear evidence for modes with angular degree / = 3. Measuring the small separation between I = 
and I = 2 allows us to plot the so-called C-D diagram of 5vq2 versus Aza The small separation 5voi of 
I = 1 from the midpoint of adjacent I = modes is negative, contrary to the Sun and solar-type stars. 
The ridge for I = 1 is notably broadened, which we attribute to mixed modes, confirming theoretical 
predictions for low-luminosity giants. Overall, the results demonstrate the tremendous potential of 
Kepler data for asteroseismology of red giants. 
Subject headings: stars: oscillations 
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1. INTRODUCTION 

Studying solar-like oscillations is a powerful way to 
probe the interiors of stars (e. g., iBrown & Gilliland 
1 19941 lChristensen-Darsgaardl2004l ) . Oscillations m mam- 
sequence stars and subgiants have been measured us- 
ing ground-based spectroscopy (for a recent review see 
lAerts et "all 120081 ) and more recentl y by space missions 
such as CoRoT (e.g., iMichel et afl 120081 ). Red giants 
oscillate at lower frequencies and so demand long and 
preferably uninterrupted time series to resolve their os- 
cillations. 

The first indications of solar-like oscillations in G 
and K giants were based on ground - based observa- 
tions in radial velocity (see [Merline! 119991 and ref- 
erenc es therein; iFrandsen et al.l l2002t IDe Ridder et al.l 
I2006Q and photometry (|Stello et a l. 2007), and on 
space-based photometry from the Hubble Space Tele 



scope (ggT;lEdmonds fc Gillilandl fl996l: iGillilandl 
Stcllo & Gillilandl 120091) . WIRE dBuzasi et al 



2008 



2000; 



Retter et all 120031 : iStello et al.l 120081). Microvariability 



and Oscillations of S tars (MOST; iBarban et al.l 120071 : 
iKallinger et al.l l2008alfbTl . and Sola r Mass Ejection Im- 
aaer (SMEI; iTarrant et all I2007D . A major break- 
through came from observations over about 150 days 
with the CoRoT space telescope, which produced clear 
detections in numerous stars of both radial and non- 
radial oscillations in t he frequency range 10-100 /xHz 
(IDe Ridder et al.ll2009t iHekker et al.l[2009t ICarTTer et al.l 
|2010[) . This increased the number of known pulsating G 
and K giants to nearly 800, and enabled the first sys- 
tematic study of oscillations i n a large population of 
red clump stars (IHekker et al.l 120091: IMiglio et all 120091 : 
IKallinger et al.l I2010I ). Note that the red clump com- 
prises core He-burning stars and is the metal-rich coun- 
terpart to the horizontal branch. 

Detecting oscillations in lower-luminosity red giants, 
which are H-shell-burning stars at the base of the giant 
branch, is very desirable for testing models of stellar evo- 
lution and also for constraining th e star-formation rate in 
the local disk (Miglio et al. 2009). However, such detec- 
tions are even more challenging than for clump stars be- 
cause the amplitudes are substantially lower. In addition, 
the oscillation frequencies are higher and are comparable 
to the orbital frequency of a satellite in low-Earth orbit. 
These factors made it difficult for CoRoT to achieve such 
detections. 

In this Letter, we present data from the aster- 
oseismology program of the NASA Kepler Mission 
(Gillil and et al.l 120101) that reveal clear oscillations in a 
large sample of red giants, extending in luminosity from 
the red clump down to the bottom of the giant branch. 

2. GLOBAL OSCILLATION PARAMETERS 

The observations were obtained over the first 34 days 
of science operations of the Kepler Mission (Ql data). 
We analyzed light curves having 29.4-minute sampling 
(long-cadence mode) for about 1500 star s listed as red 
giant s in the Kepler Input Catalog (KIC; lLatham et al.l 
I2005D . About 20% of the Ql light curves showed large 
variations (>1%) on long timescales, mostly due to M- 
giant pulsations, and are not analyzed here£3 For the 

27 One of these was subsequently found to be an oscil lating red 
giant in an eclipsing binary system (Hckkcr ct al. 2010al). 



remaining sample, we used the light curves to measure 
the global oscillation parameters, most notably the fre- 
quency of maximum power (i/ mal ) and the mean large 
frequency separation (At/). We were able to find a power 
exces s in ^1000 stars (see also Figure 5 in Gilliland ct al. 
I20l0h and measure the large separation for ^700 stars. 
In this Letter we focus on the low-luminosity red giants 
(corresponding to ^ max > 100 /iHz) for which, as dis- 
cussed in the introduction, Co RoT results have n ot yet 
been reported. As discussed bv IMiglio et al.1 (|2009), such 
stars generally have luminosities below about 30 Lq. 

Figure Q] shows the measured values of Av and ^ max 
for these stars. The symbols s how results obtaine d 
using the pipeline described by iHuber et al.l (|2009f> . 
who also described the method used to fit and sub- 
tract the background. Similar results wer e found using 



pipelin e s developed by four othe r gr oups (Hekkcr et al 
20104 iMosser & Appourchauxl 120091 : IKallinger et"al 



2010; Ma thur et alJl2010D . and the dotted lines in Fig 
ure [TJ show the region inside which ~80% of those mea- 
surements lay. A detailed comparison of the results from 
the different pipelines is ongoing. 

The tight correl ation in Figure IT] betw een Ai> and y mal 
was d iscussed by iHekker et al.1 (j2009l ) and IStello et all 
( 2009), and is seen here to extend to red giants with much 
higher ^ max . This and other correlations between global 
oscillation parameters, and their connection to funda- 
mental stellar parameters, will be discussed in future pa- 
pers. 

3. RESULTS 
3.1. Frequency analysis 

We now consider a subset of 50 of the low-luminosity 
red giants for more detailed analysis, chosen as having 
the best signal-to-noise ratios. These stars are indicated 
in Figure QJ with filled symbols. Figure [2] (left panel) 
shows power spectra for six of them, spanning the full 
range of i/ max being considered. Each power spectrum 
shows a regular series of peaks, which is the clear signa- 
ture of solar-like oscillations. 

Solar-like oscillations are p modes with high order 
and low degree, and the observed frequencies in main- 
sequence stars are reasonably well-described by the fol- 
lowing relation: 



~ Av(n + \l + e) 



l{l+l)D , 



(1) 



Here, n (the radial order) and I (the angular degree) 
are integers. The form of Equation [TJ is motivated by 
theoretical calculation s: for relatively unevolved st ars, an 
asymptotic expansion (|Tassoullll980tlGoughlll986l) shows 
that Ais is approximately the inverse of the sound travel 
time across the star, while Do is sensitive to the sound 
speed gradient near the core and e is sensitive to the 
surface layers. 

We conventionally define three small frequency separa- 
tions: 5iso2 is the spacing between adjacent modes with 
I = and I = 2, 5isi 3 is the spacing between adjacent 
modes with I = 1 and I = 3, and Sisqi is the amount by 
which / = 1 modes are offset from the midpoint between 
the / = modes on either side. If Equation [TJ holds then 
it follows that 5vq2 = 6 A), <^13 = 10-Do and Svqi = 2Do- 
We shall use Equation [TJ as a guide to the analysis of the 
observed frequencies, although its physical interpretation 
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may be open to question in the case of red giants. 

To determine oscillation frequencies for each of the 
50 stars in our subset, we needed to locate the highest 
peaks in the power spectra. We did this using conven- 
tional iterative sine-wave fitting, sometimes called pre- 
whitening or CLEAN. This involves finding the highest 
peak in the power spectrum, subtracting the correspond- 
ing sinusoidal variation from the time series and then re- 
calculating the power spectrum in an iterative process. 
We retained only those peaks that lay in the frequency 
range i^ max ± 5Av and had heights greater than 4.5 times 
the fitted background level (in amplitude). 

Damping of solar-like oscillations causes the observed 
power from each mode to be spread into multiple peaks 
under a Lorentzian envelope. If the length of the ob- 
servations is significantly greater than the typical mode 
lifetime, one therefore expects several peaks to be ex- 
tracted for each mode. In fact, the modes are only barely 
resolved (see Section [575]) and so we expect iterative sine- 
wave fitting to give good results. 

3.2. Semi-scaled Echelle Diagram 

A powerful method to study the mode frequencies of 
solar-like oscillations is to plot them modul o the large 
separ ation, in a so-called echelle diagram (|Grec et al.l 
[1983D . We have displayed the frequencies of all 50 stars 
in a single echelle diagram by using the sca ling tech- 
nique described by Beddin g fc Kieldsenl (|2010[ ) , in an ef- 
fort to carry out "ensemble asteroseismology" . The re- 
sult is shown in Figure [3]a. To make this diagram, we 
first divided the frequencies of each star by Av (that 
is, scaling them to have a large separation of unity). We 
used the echelle diagram to identify the radial modes and 
then fine-tuned the scaling factor (Az/) by up to a few 
percent to align them on a single line in the echelle dia- 
gram, shown as a sol i d line in Figure Eh. As discussed by 
iBedding fc Kieldsenl (|2010l their Method 2), this equates 
to assigning all stars the same value of e. The fact that 
all the stars can be aligned indicates this is a valid ap- 
proximation. This alignment is also clearly seen in the 
right-hand panel of Figure [2j which shows power spectra 
plotted against scaled frequencies (i.e., divided by Av). 

The ec helle diagram in Figu re [SJa differs from those 
shown by Bedding & Kieldscn (2010) in one important 
respect: only the horizontal axis has been scaled. In the 
vertical direction, we have plotted the original frequen- 
cies without scaling, which spreads out the points and 
allows us to look for variations with i/ max . Such a plot 
might be called a semi-scaled echelle diagram. 

3.3. Collapsed Echelle Diagram and Mode Lifetimes 

We can collapse the echelle diagram in the vertical di- 
rection in order to define the ridges more clearly. Fig- 
ure [3J) shows this using a simple histogram of the peaks 
in Figure [SJa, and the ridges are clearly visible. Rather 
than counting the number of peaks above a threshold, 
another approach is to sum the total power. To do this, 
we restricted each background-corrected power spectrum 
to the same frequency range as before (i/ max ±5Az/), then 
divided the frequency scale by the large separation (using 
the same fine-tuned value as above) and folded the spec- 
trum modulo unit spacing. These folded spectra were 
summed for all 50 stars and then smoothed slightly (with 



a boxcar of width O.OlAi/) to give the result shown in 
Figure [3jc. 

The clarity of the ridges in Figure [3] allows us to 
make some statements about mode lifetimes. As men- 
tioned in Section 13. 1[ damping causes each mode to 
be spread into multiple peaks under a Lorentzian en- 
velope with a full-width at half maximum (FWHM) of 
1/(7tt), where r is the mode lifetime. There has been 
considerable discussion of the lifetimes of red-giant os- 
cillations, particularly regarding whether the modes are 
sufficiently long-lived to allow their frequencies to be 



oseismologv dHoudek & Goughl 2002; 


Stcllo et all 12006: 


iBarban et al.ll2007l iDupret et al.ll2009 


). The CoRoTre- 



suits have clearly dismi s sed this worry for re d clump stars 
(|De Ridder et al.l[200l iCarrier et aLll20Toh and our Ke- 
pler data do the same for the low-luminosity stars. This 
is excellent news for asteroseismology. 

Figure [5J was constructed by aligning the radial modes 
on the vertical solid line. However, the resulting width 
of the I = ridge is not zero, which arises from at least 
three factors: (1) the resolution limit set by the duration 
of the observations (34 d), which leads to peaks in the 
power spectrum having FWHM 0.3 fiHz; (2) variations 
in Az/ with frequency in individual stars over the range 
being considered (curvature in the echelle diagram) ; and 
(3) broadening of power due to finite mode lifetimes (see 
above). The observed narrowness of the I = ridge sets 
upper limits on the second and third of these effects. 
Here, as discussed above, we are particularly interested 
in the mode lifetimes. 

The FWHM of the I = peak in Figure [3Js, found 
by subtracting the background and fitting a Lorentzian, 
is 0.033Aza This sets an upper limit on the intrinsic 
linewidths of these modes. However, so far our measure- 
ment is in terms of At/, which varies over the sample. 
The median Av is 11.5 //Hz, in which case the measured 
width equates to a lower limit on the mode lifetimes of 
about 10 days. To investigate further, we divided the 
sample into four bins according to Av and formed the 
collapsed echelle diagram for each subset, as before. We 
found the widths of the I = ridges in all four cases to be 
about the same in absolute terms (~0.4/xHz), which is 
only slightly greater than the FWHM set by the resolu- 
tion of the data (see above) . We conclude that the mode 
lifetimes in these red giants is at least 10 days, possibly 
much greater. 

3.4. Small Separations and the C-D Diagram 

The I = 2 ridge is closely parallel to the I = ridge, 
implying that 8vq2 is a nearly constant fraction of Av 
for the stars in our sample, with 8v^2 ~ 0.125Az/ (Fig- 
ure [5Jc). This frequency scaling indicates that, to a good 
approximation, the stars are homologous. 

We were able to measure the mean small separation 
8vq2 for 38 stars in our sample. Figure Bh shows the re- 
sults in a so-called C-D diagram (|Christensen-Dalsgaardl 
1988), which plots 5v 02 versus Av. The dashed line in 
Figure shows a linear fit, which has parameters 

5v m = (0.122 ± 0.006)Az/ + (0.05 ± 0.08) //Hz. (2) 

The spread of points about this relation presum- 
ably reflects a spread in stellar masses, and is better 
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seen by plotting the rati o of the two separations, as 
shown in Figure [4fr (see iRoxburgh fc Vorontsovl 120031 : 
lOti Floranes et al1T200 5; MazumdarJl2Q05() . 

The dotted line in Figure |3] marks the midpoint be- 
tween I = modes. Interestingly, the I = 1 ridge is 
centered slightly to the right of this line. This indicates 
that the average small separation Svqi in these red gi- 
ants is negative, whereas it is positive in the Sun and 
the handful of other main- sequence stars for w hich it has 
been measured: a Ce n A (iBedding et al.ll2004fl . a Cen B 
(iKieldsen et al.l 12005ft 7 P av (jMosser et al.l I2008T ). and 
t Cet (jTeixeira et al.lT2009h . A negative value for 8vq\ 
also indicates that Equation [TJ is not satisfied in detail 
for these red gia nts, a point a l ready noted for the G8 
giant HR 7349 bv lCarrier etaLl (|20Tol) . based on C0R0T 
observations. 

3.5. Mixed Modes with I = 1 

The I = 1 ridge in Figure [3] is notably broader than 
the others. We do not attribute this to shorter mode life- 
times, but rather to the presence of mixed modes whose 
frequencies are shifted by avoided crossings. Mixed 
modes occur in evolved stars and have the characteristics 
of p modes in the envelope of the star an d of g modes in 
the interior (e.g.. lAizenman et al.lH977l ). A typical ex- 
ample is the star KIC 5356201 (see Figure [5] and hlled 
blue symbols in Figure [3]a) . The I = and 2 modes are 
quite close to their respective ridges, whereas there is 
much more scatter about I = 1. A more extreme exam- 
ple is the star KIC 4350501 (see Figure [2 and filled red 
symbols in Figure [3]a) , which shows a particularly high 
peak at 87/iHz, well below v max . The proximity of this 
peak to the I = 1 ridge indicates that it is likely to be a 
mixed mode. A similar situation is seen in the F5 sub- 
giant Procyon, for which ground-based velocity observa- 
tions show a narrow peak well below v max that lies close 
to th e I = 1 ridge in the echelle diagram (Bedding et al. 
2010). The large peak height can be explained by the fact 
that mixed modes are expected to have longer lifetimes 
(smaller linewidths) than pure p modes because they 
have larger mode inertias (e.g.. IChristensen-Dalsgaardl 
12001 . 

Another good example is KIC 5006817, shown in Fig- 
ure [SJ where we clearly see multiple I — 1 peaks in each 
order (see also KIC 9904059 = KOI 145 in Figure 1 of 
IGilliland et al.l I2010T) . Each of these clusters of peaks 
is reminiscent of a single mode that is heavily damped 
(see Section |3.3[) . However, as mentioned above, mixed 
modes are expected to have longer lifetimes than pure p 
modes and hence to produce narrower - not broader - 
peaks. Indeed, this phenomenon of multiple I = 1 peaks 
is probably what led previous authors to report short 
mode lifetimes in red gi ants (see Sectionl3.3l) . Guided by 
the theoretical work of IDupret et al.l (|2009l ). we instead 
interpret the observed properties of the / = 1 ridge in 
t erms of mixed modes . 

IDupret et al.l ()2009f ) discussed theoretical amplitudes 
of solar-like oscillations for a range of models on the 
red giant branch. Their models show mostly regular 
frequency patterns for evolved giants (f ma x < 80^tHz), 
but predict more complex spectra for low-luminosity gi- 
ants. This is mainly explained by the less efficient ra- 
diative damping in the cores of these stars, leading to 
a stronger interaction of p-mode and g-mode cavities, 



and by the less efficient trapping due to the smaller 
density contrast between the core and the envelope. 
As a consequence, more non-radial modes could be ob- 
served than only those trapped in the envelope. This 
is particularly expected to affect I = 1 modes because 
the evanescent region just above t he H-burning shell 
is smaller for these modes (see also Dzicmbowsk i et all 
1200 It IChristensen-Dalsgaardl 120041 ). Our observation 
of a broadened I — 1 ridge with a small cluster of 
modes per order sup ports the theoretical predictions by 
IDupret et al.l (|2009l) . Observations of frequency spectra 
in a sample of higher luminosity red giants would be a 
desirable follow-up study. 

3.6. Detection of I = 3 Modes 

In addition to ridges corresponding to I = 0, 1 and 2, 
we see clear evidence in Figure [3] for modes with 1 = 3. 
The detection of these modes with photometry is difficult 
because of the very low amplitudes that result from ge- 
ometric cancellation. The only previous report of solar- 
like oscillations with I = 3 from photometry (except for 
the Sun) is a probable detection using C0R0T in the G O 
mam-sequence star HD 49385 (jDeheuvels et al.l l20Tol ). 
The ability to detect I = 3 modes from Kepler obser- 
vations of red giants is a significant bonus. 

In Figures [3)3 and b there are 12 peaks that fall along 
the ridge that we identify with I = 3, and these arise 
from 10 stars in the sample of 50. These are presumably 
peaks that, due to the stochastic nature of the excitation, 
happen to rise above our detection threshold in these 
observations. A more objective approach is to sum the 
power from all stars, as shown in Figure \5jp, and we again 
see clear evidence for the I = 3 ridge. 

How does the position of the I = 3 ridge compare with 
expectations? As mentioned in Scction f3.il it is conven- 
tional to measure I = 3 relative to / = 1, via the small 
separation SV13. However, given that the position of I = 1 
for these red giants is anomalous, it seems sensible to in- 
stead measure I = 3 relative to / = 0. We therefore 
suggest using a new small separation, Svqz, defined (by 
analogy to Si^oi ) as the amount by which the I = 3 modes 
are offset from the midpoint between the / = modes on 
either side. If follows that 5vo3 = Svoi + <5z^i3 . From Fig- 
ure^ we see that SV03 ~ 0.28Aza If the asymptotic rela- 
tion held (Equation[lJ , we would expect S1S03 = 12Dq. In 
that case the ratio between 61^03 and Svq2 would be 2.0, 
whereas we observe a ratio of 2.2. These results should 
provide valuable tests for theoretical models. 

4. CONCLUSIONS 

These results, based on the analysis of only 34 days 
of data, show the tremendous potential of Kepler long- 
cadence data for the study of solar-like oscillations in 
red giants. As more data become available, the longer 
time series will allow detailed studies of stars covering 
the whole range of evolutionary states along the red gi- 
ant branch and provide stringent observational tests for 
theories of stellar structure and evolution in this part of 
the H-R diagram. 

We gratefully acknowledge the entire Kepler team, 
whose outstanding efforts have made these results pos- 
sible. Funding for this Discovery mission is provided 
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Fig. 1. — Large frequency separation versus frequency of maximum power for 78 low-luminosity red giants. Filled symbols indicate the 
50 stars that were selected for more detailed analysis. The dotted lines enclose the region containing ~80% of measurements by the five 
pipelines (see text). 
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Fig. 2. — Left: Pow er s pectra of six representative low-luminosity red giants. Right: The same power spectra plotted against scaled 
frequency (see Section |3.2|I . The dotted lines are equally spaced, having u nit separation and b eing aligned with the I = modes. Stars are 
labelled with identification numbers from the Kepler Input Catalog (KIC; Latham ct al. 2005). 
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Fig. 3. — (a) Semi-scaled echelle diagram for the sample of 50 low-luminosity red giants (see Section |3. 2D . Symbols indicate extracted 
frequencies, with symbol sizes proportional to am plitu de. Filled blue and red symbols are frequencies of the stars KIC 5356201 and 
KIC 4350501, respectively (see Figurc[2]and Scction l3.5|l . The value of Au for each star was fine-tuned in order to align the £ = modes on 
the left-most vertical dotted line. The other dotted line lies exactly 0.5 away, marking the midpoint between £ = modes. (f>) Histogram 
of the poin ts shown in panel a. (c) Folded and scaled power spectra, collapsed over all stars in the sample and smoothed slightly (see 
Section l3.2ll . The arrows show the definitions of the s mall s eparations Svqi, Suo2, an d Svi3, with the sign convention that left-pointing 

arrows correspond to positive separations (see Section l3.lt . 
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Fig. 4. — (a) The so-called C-D diagram ( Christcnscn-DalsKaard 1988) of Suo2 versus Av for the 38 stars in our sample for which both 
quantities were reliably determined. The dashed line shows a linear fit. (6) The ratio &UQ2/ versus Av for the same stars, to better show 
the scatter about the trend. 
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Fig. 5. — Power spectrum and echelle diagram for the star KIC 5006817, illustrating multiple 1 = 1 peaks per order. Note that the highest 
peak in the power spectrum extends beyond the plot limits (to 8400 ppm 2 ). 
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